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ABSTRACT 

Identifying liquid water on the surface of planets is a high priority, as this traditionally defines 
habitability. One proposed signature of oceans is specular reflection ("glint"), which increases the 
apparent albedo of a planet at crescent phases. We post-process a global climate model of an Earth- 
like planet to simulate reflected lightcurves. Significantly, we obtain glint-like phase variations even 
though we do not include specular reflection in our model. This false positive is the product of two 
generic properties: 1) for modest obliquities, a planet's poles receive less orbit-averaged stellar flux 
than its equator, so the poles are more likely to be covered in highly reflective snow and ice, and 2) we 
show that reflected light from a modest-obliquity planet at crescent phases probes higher latitudes than 
at gibbous phases, therefore a planet's apparent albedo will naturally increase at crescent phase. We 
suggest that this "latitude-albedo effect" will operate even for large obliquities: in that case the equator 
receives less orbit-averaged flux than the poles, and the equator is preferentially sampled at crescent 
phase. Using rotational and orbital color variations to map the surfaces of directly imaged planets 
and estimate their obliquity will therefore be a necessary pre-condition for properly interpreting their 
reflected phase variations. The latitude-albedo effect is a particularly convincing glint false positive for 
zero-obliquity planets, and such worlds are not amenable to latitudinal mapping. This effect severely 
limits the utility of specular reflection for detecting oceans on exoplanets. 



1. INTRODUCTION 

The traditional habita ble zone (HZ) is defi ned in terms 
of surface liquid water (|Kasting et al.|[l993() . Three dis- 
tinct methods have been proposed to search for liquids 
on the surface of a planet: 

Rotational color variability (D>rd et al.l l200l|: 
Cowan et al.ll2009t iKawahara fc Fuiiill2010l : lCowan et al.l 
20111 : iKawaha ra fc Fuiiil 1201 1[ ): oceans are darker and 



have different colors than other surface types on Earth, 
so the time variations in color of a spatially unre- 
solved planet can betray the presence of liquid water 
oceans. This method relies on there being longitudinal 
inhomogeneities in the planet ' s surf a ce com position. 

Polarization (jZugger et al.l 120101 1201 ID : oceans are 
smoother than other surface types and thus polarize 
light. For idealized scenarios, the phase variations in 
polarization are significant, but the same authors found 
that in practice the effect of oceans is masked by Rayleigh 
scattering, clouds and aerosols. Observations of polar- 
ized Earthshine, however, imply that rotational varia- 
tions in polarizatio n may be useful in detecting oceans 
(jSterzik et alJl20"I3 ). 

Specular reflectio n (jWilliams fc Gaidosl 120081 : 
iRobinson etld1l2l)Tol ): oceans are also able to specu- 
larly reflect light, especially at crescent phases. The 
signal-to-noise requirements for phase variations are 
not as stringent as for rotational variations since the 
integration times c an be much longer: w eeks instead 
of hours. However, IRobinson et al.1 (|2010l ) showed that 
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clouds not only mask underlying surfaces, but forward 
scattering by clouds mimics the glint signal at crescent 
phases, while atmospheric absorption and Rayleigh 
scattering mask the glint signature. They proposed 
using near infrared opacity windows to search for glint, 
but this would only be possible if the effects of clouds 
could be accurately modeled for exoplanets. 

Clearly, any method attempting to probe the planetary 
surface will be impeded by clouds, but the three tech- 
niques above have been shown to work in limited empiri- 
cal and/or numerical experiments for planets like Earth, 
which is roughly half covered by clouds. Such tests are 
necessary but not sufficient. While previous authors have 
noted that seasonal changes in snow, ice and cloud cov- 
erage le ad to 15-20% modulations in a planet's apparent 
albedo (IWilliams fc Gaidosl [20081 IQaklev fc Cashll200l 
IRobinson et al.ll2010f) . they under-appreciated the im- 
portance of viewing geometry. 

In this Letter, we focus on the specular reflection 
method and show that the generic tendency of HZ plan- 
ets to have snow and ice in their least-illuminated regions 
naturally leads to a false positive for ocean glintO This 
is because the viewing geometry changes throughout an 
orbit such that the least-illuminated latitudes are pref- 
erentially sampled at crescent phases. We call this the 
latitude-albedo effect. 

2. METHODS 

2.1. Global Climate Model 

We calculate reflected phase variations for a vari- 
ety _of_viBwing geome tries using a simulation presented 
in iVoigt et "all (|2011[ ) generated by the global climate 
model (GCM) ECHAM5 /MPI-OM. It simulates the cir- 

5 We demonstrate the latitude-albedo effect using a model with 
an ocea n, but nearly dry planets develop similar latitudinal albedo 
profiles ijAbe et al. 
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culation of the atmosphere and ocean and has inter- 
active sea-ice and c louds; a det a iled d escription of the 
GCM is given in iVoigt et al] (|2011l ) and references 
therein. ECHAM5/MPI-OM includes a shortwave radia- 
tion model that calculates the top-of atmosphere incident 
(FjJ and reflected (Ff) shortwave radiative flux at every 
time step and grid location of the model. 

The model has a 3.75° horizontal resolution and does 
not allow radiative transfer between columns. The short- 
wave scheme includes absorption and scattering. For- 
ward scattering from aerosols and cloud particles is in- 
corporated by asymmetry factors and using the delta- 
Eddington approximation. The model accounts for the 
increased atmospheric path length at grazing angles of 
incidence, but surface reflection does not depend on so- 
lar zenith angle and there is no specular reflection in our 
GCM. 

This simulation uses a modern-Earth solar constant 
of 1367Wm~ 2 , pre-industrial atmospheric greenhouse 
modern-Earth orbit (23.5° obliquity, 1.7% ec- 
centricity) and continent configuration meant to repre- 
sent Earth 635 Ma before present. The simulation is 
initiated from a climate state with polar ice caps and 
remains in such a state throughout. It was run to equi- 
librium, then for an additional ten years, with shortwave 
radiation outputs saved every two hours. The planet has 
a mean surface albedo of 0.244 and is 57% covered in liq- 
uid water, making it habitable by definition. Cloud cover 
is 0.654 and the time-averaged planetary albedo is 0.351, 
both slightly higher than for simulations of pre-industrial 
modern Earth. 



We first consider the bolometric spherical albedo of the 
planet, or Bond albedo, 



A B (t) = 



' F t (e,(f),t)dn 



Ft 



(1) 



where § dVL is the surface integral over the entire sphere. 

The Bond albedo exhibits seasonal variations of 6% 
(Figure [2J . These are largely driven by the meridional 
gradient in albedo (right panel of Figure[T]) and are much 
smaller that the latitude-albedo effect. 

2.3. Disk-Integrated Reflected Light 

We now calculate light curves as they would appear to 
a distant observer, assuming diffuse (a.k.a. Lambertian) 
reflection. We completely neglect any specular reflection 
that could produce a glint spot, yet, as we show below, 
we still find that the apparent albedo exhibits the phase 
variations one would expect for glint. 

We define the normalized weight as W(9, 4>, t) = 
VFi/ § VF^dQ, where V(9, <j>, t) is the visibility of a given 
region of the planet for a given observer (V equals one at 
the sub-observer point, drops as the cosine of the angle 
from the sub-observer point, and is zero on the far side of 
the planet). For a given viewing geometry, W quantifies 
the sampling of different regions on the planet. 

The apparent albedo is the disk-averaged albedo of the 
planet, weighted by illumination and visibility: 



A*(t) 
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Fig. 1. — Ten-year averaged map of planetary albedo for the sim- 
ulation. Green lines indicate the coast of the equatorial continents; 
blue lines show the annually averaged extent of sea-ice. The right 
panel shows the zonally (longitudinally) averaged albedo profile. 

Since we use top-of- atmosphere fluxes, all of the albe- 
dos in the remainder of this Letter are planetary albedos, 
as opposed to surface albedos. We define A(9, (f>, t) = 
F^-(9,<p,t)/Fi(9,(t),t), where 9 and 4> are planetary lat- 
itude and longitude, respectively. In Figure [T] we plot 
A = F^/F^, were overbars indicate time-averages. 

We compute lightcurves with 2 hr time resolution, in- 
cluding rotational variations, as well as changes in cloud 
cover on diurnal, seasonal and inter-annual timescales. 
For the purposes of this Letter we consider 24-hr inte- 
gration times for the observations, which average over 
the rotational variations. 

2.2. Bond Albedo 



3. RESULTS 

In the top panels of Figure [2] we show the phase vari- 
ations in apparent albedo for two generic observing ge- 
ometries (solid black lines). We gray-out the phases for 
which the planet would be inside the inner working an- 
gle (IWA) of a high-contrast imaging mission and there- 
fore unobservable. The IWA depends on the details of 
a mission's design and is a function of wavelength and 
distance. We gray-out phases within 45° of full or new 
phase, corresponding to an IWA of 71 mas for a system 
at 10 pc. 

The phase variations in apparent albedo are shown 
with solid black lines in the top panels of Figure [2J De- 
spite the fact that our model includes only diffuse reflec- 
tion, the variations have approximately the same shape 
and amplitud e as glint for an Earth-lik e geography and 
cloud cover ({Williams fc Gaidosl 120081 : iRobinson et al.1 

|2oToT) . 

For the left geometry of Figure [2J the planet / star flux 
ratio is 2.4 x 10 -10 and 0.6 x 10~ 10 at superior and in- 
ferior conjunction, respectively. For an Earth-twin at 
10 pc, a 10 m telescope with 5% throughput could ob- 
tain 1% photometry in a 100 nm optical band in 6 days 
of integration at i nferior conjunction, w hen the planet is 
faintest (following iRoberge et al.ll2012f ). 

3.1. Dominant Latitude 

The phase variations can be understood in terms of 
the dominant latitude and the planet's latitudinal albedo 
profile. The dominant latitude is the latitude most sam- 
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Fig. 2. — Top: Phase variations in apparent albedo (solid black) compared to Bond albedo (blue) for two viewing geometries. The dashed 
black lines show the apparent albedo for a static, time-averaged map. The left and right edges of the plot correspond to crescent phases. 
The grey regions in the right panel denote the phases that will not be accessible because the planet will be inside the inner working angle 
(IWA) of the telescope; the geometry shown at left is sufficiently face-on that the planet is never within the IWA. Bottom: Phase variations 
in the dominant latitude, 0dom> f° r the same two viewing geometries. The orbital inclination, i, is the angle between the planet's orbital 
angular momentum and the line of sight (i = 0° for a face-on orbit, i = 90° is edge-on); the sub-observer latitude is denoted by o b s - 
The time of periastron, t per i, determines the orientation of the (slightly elliptical) orbit with respect to the line of sight. Northern winter 
solstice nearly coincides with periastron, as on present-day Earth. 



pled for a given viewing geometry ([Cowan et al.ll20lil) : 



(0) = Iw{9,(f>,t)9dn. 



(3) 



In the bottom panels of Figure we show the sub- 
stellar, sub-observer, and dominant latitudes for the 
same two geometries used in the top panels. Significantly, 
the dominant latitude may be poleward of both the sub- 
stellar and sub-observer latitudes, but this only occurs 
at crescent phases (near inferior conjunction). 

The latitudinal albedo profile is shown in the right 
panel of Figure [TJ Albedo is lowest in the tropics, so 
the planet's apparent albedo will tend to be greater than 
its Bond albedo, because the dominant latitude is typi- 
cally poleward of the sub-stellar latitude. Furthermore, 
variations in the dominant latitude over the course of 
a year couple with the latitude dependence of albedo to 
generate phase variations in the planet's apparent albedo 
(the latitude-albedo effect). 

We estimate the magnitude of the latitude-albedo ef- 
fect by replacing the actual time-variable albedo of the 
planet with its ten- year average (Figure [1]) when calcu- 



lating apparent albedo (Eq. [2]). Since the 10-year albedo 
map averages over diurnal and seasonal variations in 
albedo, this exercise isolates the latitude-albedo effect. 
The approximate albedo variations are shown by the 
dashed black lines in the top panels of Figure [2j 

Clearly the latitude dependance of albedo is the dom- 
inant factor enhancing the apparent albedo at crescent 
phase in our simulation. The discrepancies between the 
solid and dashed black lines are due to second-order ef- 
fects. For example, seasonal ice and snow increase the 
apparent albedo during Boreal winter and spring (top- 
left panel of Figure [2]), while the dependence of albedo 
on solar zenith angle provides an additional enhancement 
of apparent albedo at crescent phases (top-right panel 
of Figure [2]). At large zenith angles, GCM pixels with 
clouds overlying dark surfaces (water or land) have en- 
hanced albedos because the longer path-length through 
the atmosphere leads to more scattering. This increases 
the brightness of the day-night terminator at low lati- 
tudes. Diurnal cycles in cloud coverage could also cause 
reflected p hase variations (f or Earth these patterns are 
weak; lHartmann et al.lll991l ) . 
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3.2. Statistical Analysis of Geometry 

The latitude-albedo effect is a confounding factor for 
the photometric detection of specular reflection. To be a 
convincing false-positive, however, the albedo must not 
only increase at crescent phases (true of both geometries 
in Figure [2]), but must also be symmetrical at waxing and 
waning phases (true for the left geometry, but not the 
right). The critical variable affecting these observables is 
the orientation of the planet's axial tilt with respect to 
the observer: the obliquity phase. 

We perform a statistical analysis of inclinations, obliq- 
uities, and obliquity phases to determine which viewing 
geometries might be convincing false positives. We can- 
not predict how snow and ice coverage will vary as a 
function of obliquity without rerunning computationally 
expensive climate simulations, but we track dominant 
latitude as a function of orbital position for each geome- 
try. In particular, we consider (\8\) = § W\d\dVL since we 
only care about distance from the equator, not whether 
a region is in the Northern or Southern hemisphere. 

The most important quantities are 1) the amplitude 
of the excursions in dominant latitude: at inferior 

conjunction minus that at superior conjunction (this un- 
derestimates the amplitude for asymmetric excursions), 
and 2) the asymmetry in (\9\) about superior conjunc- 
tion, which we estimate as at waning quarter minus 
at waxing quarter. In Figure [3] we adopt the most 
likely orbital inclination, 60°, and plot the amplitude of 
dominant latitude excursions versus the asymmetry of 
those excursions. 




Fig. 3. — Excursions in dominant latitude, (\Q\), for planets with 
60° orbital inclination. Each color shows the locus for a single 
obliquity, ©, and 1000 different obliquity phases; the black as- 
terisk shows the zero-obliquity case. Positive amplitudes (top of 
plot) indicate that the regions probed at crescent phase are pole- 
ward of those sampled at gibbous phase; this is always the case 
for planets with modest obliquities, © = 0° and 30°. Negative 
amplitudes (bottom of plot) indicate that the regions sampled at 
crescent phases are equatorward of those probed at gibbous phases; 
this is always the case for © = 90° . Non-zero asymmetry (left and 
right edges of plot) indicate that the latitudes sampled at waning 
phases are different from those sampled at waxing phases. The 
white regions correspond to (|9|)-variations that are liable to be 
convincing glint false positives: the change at crescent phases will 
be greater than the asymmetry between waxing and waning phases. 
The fraction of planets falling in this "dangerously convincing" area 
is 100%, 59%, 30% and 29% for obliquities of 0°, 30°, 60°, and 90°, 
respectively. 



Recall that for obliquities less than 53.9° the equator 
receives more orbit-averaged insolation than the poles, 
while at greater obliquities the poles receive more flux 
than the equator, albeit with large seasonal variations. 
In nearly all cases shown in Figure |3l the regions probed 
at crescent phase receive less orbit-averaged insolation 
than those probed at gibbous phases, which will lead 
to larger apparent albedos at crescent phase if these re- 
gions have year-round snow. The exception is = 60°, 
for which the latitudes sampled at crescent phase may 
be poleward or equatorward of those probed at gibbous 
phases; but a planet with such an obliquity has relatively 
uniform orbit-averaged insolation as a function of lati- 
tud e and may not have year-round "cold" regions (but 
see lAbe et alJ[2TjTTI) . 

The glint false positive is most extreme and pernicious 
for planets with zero obliquity because the dominant lat- 
itude at crescent phases is nearly 40° poleward of that 
at gibbous phases, and the excursions are perfectly sym- 
metrical about superior conjunction (black asterisk at 
top of Figure [3]) . For planets with non-zero obliquity, 
the latitude-albedo effect best mimics the glint signature 
(no asymmetry) when the orbital/viewing geometry is 
such that solstices coincide with conjunctions. 

The same effect occurs for edge-on inclinations (not 
shown), but in that case the dominant latitude only 
increases significantly at extreme crescent phases (e.g., 
within ~ 6° of inferior conjunction, for i = 85°). 

4. DISCUSSION AND CONCLUSIONS 

The latitude-albedo effect operate if the regions of the 
planet receiving the least orbit-averaged flux have the 
greatest albedo due to inhomogeneous surface or cloud 
c over. 

iWilliams fc Pollard! (|2003| ) simulated the climate of 
high-obliquity Earth analogs and did not find year-round 
equatorial snow, but they speculated that such cover- 
age would occur for lower CO 2 levels, or higher-elevation 
continents. lAbe et all ([201 1[ ) found that planets with no 
large surface liquid water reservoirs store their condensi- 
bles at their coldest latitudes, regardless of obliquity. 

Even if no surface snow/ice exists on a planet, or if the 
surface is obscured by a thick atmosphere, the latitude- 
albedo effect will still operate if the coldest regions of a 
planet are most cloudy. This is the case in our simula- 
tion, but may not be in general. Meridional gradients 
in cloudiness could be important for planets with surface 
temperatures too hot for snow and ice, due to low albedo, 
high insolation, high concentrations of greenhouse gases, 
or thicker atmospheres — in short, planets at the inner 
edge of the HZ. 

If a planet's apparent albedo increases at crescent 
phases, there are two fundamentally different explana- 
tions: 1) the scattering phase function of the planet is for- 
ward peaked due to, for example, atmospheric Rayleigh 
scattering, Mie scattering from clouds, or specular reflec- 
tion from water, or 2) there are more scatterers in the 
scene, either intrinsically, as with the latitude-albedo ef- 
fect, or because of path length effects. 

If a planet's apparent albedo exceeds unity, then — by 
definition it must be scattering non-diffusively, but this 
only occurs at or with in the IWA for Earth-like models 
(|Robinson et al.ll2010f) . More importantly, the albedo- 
radius degeneracy for directly imaged planets makes it 
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impossible to put albedo on an absolute scale. 

Glint polarization fraction peaks at quadrature and 
cannot directly distinguish between sources of brighten- 
ing at crescent phase, although polarimetry could poten- 
tially detec t surface water at other phases (IZugger et al.l 
l20la.l20Tl. 

iRob inson et a l. (2010) noted that glint measurements 
should be performed at near-infrared (NIR) wavelengths 
in order to avoid Rayleigh scattering, and in opacity win- 
dows in order to avoid atmospheric absorption. The 
reflectance of the glint spot is essentially the same at 
NIR wavebands as in the optical, while snow becomes 
less reflective at these longer wavelengths. The optical- 
NIR color of snow depends sensitively on grain size, with 
larger grains being redder ()Warrenl [19821 . It is therefore 
possible that measurements of reflected phase variations 
at both optical and near-infrared wavelengths could help 
disentangle the effects of ocean glint from that of large- 
grained snow, provided that Rayleigh scattering does not 
obscure the surface at the shorter wavelength. 

To summarize, detecting ocean glint for an Earth-like 
planet based solely on multiband phase variations is only 
possible if the effects of clouds and snow can be properly 
modeled for extrasolar planets. For the remainder of this 
Letter we assume that the snow/glint degeneracy will not 
be resolved, and examine how the latitude-albedo effect 
is likely to affect exoplanet characterization. 

High contrast imaging missions will be able to mon- 
itor the apparent albedo of terrestrial planets in the 
HZ of nearby Sun-like stars. In the nea r-term, pairing 
the James Webb Space Tel escope (JWST;lGardner et al.l 
2006) with a starshade (jSoummer et al.l 12010ft might 
achieve the same goal. 

If a directly imaged terrestrial planet is observed for 
only half an orbit and its apparent albedo is greater at 
crescent than at gibbous phase, no conclusion can be 
drawn. If the observed albedo variations are monitored 
for a full orbit and are symmetrical, it will mean a) the 
planet has zero obliquity, b) the planet's solstices coin- 
cide with conjunctions, or c) the latitude-albedo effect 
is unimportant. If, on the other hand, the lightcurve is 
strongly asymmetrical, one will suspect that the latitude- 
albedo effect is at play, but will be unable to correct for 
the effect, because of the unknown obliquity and obliq- 



uity phase of the planet. Note that asymmetry in appar- 
ent albedo may also be due to seasonal changes in the 
planet's intrinsic Bond albedo, but we found this to be a 
minor effect in our simulation. 

Fortunately, high signal-to-noise, high-cadence re- 
flected light measurements of an imaged exoplanet can 
yield precisely the variables one needs to correct for 
the latitude-albedo effect. This involves: 1. deter- 
mining the rotation rate of t he planet (tested in sim- 
ulations bv IPalle et all [20081: lOaklev fc Gashl l2009h . 2. 
making rotational albedo maps of th e planet (tested us- 
ing EPOXI Ear th observations by iCowan et al.1 120091 : 
iFuiii et aT1l2011l ). and 3. doing so at a variety of orbital 
posit i ons (tested in simulations b y iKawahara k, FuiiH 
l20loll20ll[Fuiii fc Kawaharal[20Tl . 

It is therefore possible to properly interpret re- 
flected light phase variations of imaged planets at cres- 
cent phases provided that one obtains high signal-to- 
noise, high-cadence measurements at gibbous phases 
(for signal-to-noise consi derations, see lCowan et all2009t 
IKawahara k Fuiii|[2010h . 

Zero-obliquity planets are the worst case for three rea- 
sons: 1) variations in dominant latitude are large and 
symmetrical, 2) poles receive the least orbit-averaged 
flux, and are therefore most likely to harbor year-round 
snow, and 3) these planets are not amenable to latitu- 
dinal mapping, because the sub-stellar point is always 
equatorial. Although planets with negligible obliqui- 
ties will likely be a minority of directly-imaged worlds, 
they are the norm in the HZ of low-mass stars. The 
latitude-albedo effect will therefore be an important glint 
false positive for missions characterizing the reflected 
phase variations of temperate planets orbiting nearby M- 
Dwarfs. 
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